Magnetic-Field-Induced Stabilization of Non-Equilibrium Superconductivity 



J. T. Peltonen, 1 J. T. Muhonen, 1 ' 2 M. Meschke, 1 N. B. Kopnin, 1 ' 3 and J. P. Pekola 1 

1 Low Temperature Laboratory, Aalto University, P.O. Box 13500, FI-00076 A ALTO, Finland 
2 Department of Physics, University of Warwick, CV4 7AL, UK 
3 L. D. Landau Institute for Theoretical Physics, 117940 Moscow, Russia 
(Dated: August 9, 2011) 

A small magnetic field is found to enhance relaxation processes in a superconductor thus stabilizing 
superconductivity in non-equilibrium conditions. In a normal- metal (N) - insulator - superconductor 
(S) tunnel junction, applying a field of the order of 100 fiT leads to significantly improved cooling of 
the N island by quasiparticle (QP) tunneling. These findings are attributed to faster QP relaxation 
within the S electrodes as a result of enhanced QP drain through regions with locally suppressed 
energy gap due to magnetic vortices in the S leads at some distance from the junction. 



PACS numbers: 74.50. +r, 74.25.Ha, 73.40.Rw 

In this Letter, we report an observation that appears 
counterintuitive at first: a small magnetic field is found to 
stabilize superconductivity under quasiparticle (QP) in- 
jection. In our experiment, the cooling power of normal- 
metal (N) - insulator (I) - superconductor (S) tunnel 
structures is enhanced in perpendicular magnetic fields 
B± ~ 100 /iT = 1 G. Measured maximum tem- 
perature drop ST relative to a starting bath tempera- 
ture To = 285 mK exhibiting this behavior is shown in 
Fig. [I] (a). The improvement is unexpected, as in general 
the effect of a magnetic field is to suppress superconduc- 
tivity. Electronic cooling in NIS junctions in the presence 
of magnetic fields in both perpendicular and parallel ori- 
entations has been studied before [1 , but only in higher 
fields where the cooling power was already reduced due 
to diminishing superconducting energy gap A. On the 
other hand, the creation of magnetic vortices [2] has been 
shown to enhance QP relaxation in superconducting alu- 
minum, as the QPs become trapped and thermalize in the 
regions of reduced A [3 . Here we demonstrate that the 
additional relaxation channel due to enhanced QP drain 
through regions occupied by magnetic vortices enhances 
the superconducting performance of S leads and improves 
the electronic cooling in NIS junctions. This can be of rel- 
evance in superconducting qubits [IHZ] , resonators [8] and 
in hybrid SINIS turnstiles [9] in reducing the effects from 
nonequilibrium and residual QPs arising due to drive and 
microwave radiation from the environment. Moreover, 
improved relaxation caused by vortex creation in the S 
leads can partially explain the "re-entrant superconduc- 
tivity" observed in Zn and Al nanowires [TO] [11] . In the 
present case, as sketched in Fig.[l](a), vortex formation in 
the S electrodes away from the NIS junction improves re- 
laxation of the injected QPs and leads to enhanced cool- 
ing of the N island. In higher fields, vortices move closer 
to the junction, deteriorating the cooling power. 

In an NIS junction, A acts as an energy filter for the 
tunneling QP [T2HT5] . At low temperatures k^T <C A 
and for bias voltages eV < A across the junction, the 
electrons in the N electrode cool considerably below the 




FIG. 1. (color online) (a) Maximum temperature drop ST 
in an optimally biased SINIS cooler in a perpendicular mag- 
netic field B±, at the bath temperature To = 285 mK. The 
sketches show the S electrode geometry and qualitative vortex 
configurations at B± < 2 G and at a value of B± beyond the 
optimum point. The area inside the green dashed rectangle 
corresponds to that in the micrograph below, (b) Scanning 
electron micrograph of a typical structure, together with the 
measurement scheme (see text for details) . A Cu island (red) 
is contacted to four superconducting Al electrodes (blue) via 
Al oxide tunnel barriers for thermometry and temperature 
control. Replicas of each structure are visible due to the fabri- 
cation involving two-angle shadow evaporation of the metals. 



phonon temperature by hot QP extraction. The effect 
can be made more pronounced in a symmetric double 
junction SINIS structure with a small N island contacted 
to S leads via two NIS junctions [15 , allowing to con- 
struct practical solid-state refrigerators for cooling thin- 
film detectors to temperatures close to 100 mK [16] [17] . 
The performance of actual devices depends crucially on 
the relaxation of the QPs that are injected into the S elec- 
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trode, as the superconductor overheating diminishes the 
cooling power at an NIS junction because of enhanced QP 
backtunneling. The excess QP density close to the junc- 
tion can be diminished by fabricating the S electrodes 
very thick [16], or covering them partially by a layer 
of normal metal that acts as a QP trap [T8H20] . The 
QP population is typically modeled in terms of a diffu- 
sion equation, describing their recombination retarded by 
phonon retrapping, and other loss mechanisms [2TH25] . 
Converting the excess density into an effective, position- 
dependent temperature T(x) [26 , 27 , one finds that at 
phonon temperatures k^T <C A the S leads can be over- 
heated on a length scale ranging from tens of micrometers 
to a millimeter, as the electron-phonon relaxation and 
electronic heat conduction are exponentially suppressed 
compared to their normal state values [22 } [28 ] . 

Here we present data from one of several measured 
symmetric SINIS structures similar to that shown in 
Fig. [l](b), fabricated at different times and electrically 
characterized in a dilution refrigerator down to 50 mK 
bath temperature. The same qualitative behavior was 
observed in all structures with the same geometry. A cop- 
per island of area ^4n — 2.7 x 0.7 fim 2 is contacted by four 
over lap- type Al/Al-oxide/Cu NIS junctions. The Al elec- 
trodes with zero-temperature energy gap Ao ^ 210 fieV 
become superconducting below Tq c± 1.4 K. Compared 
to the two small (probe) junctions in the middle, the two 
outer (cooler) junctions at the ends of the island have 
larger overlap area and therefore lower normal state tun- 
nel resistance Rt — 1.1 kf2 each. The structures were 
fabricated on an oxidized silicon substrate by standard 
electron beam lithography and two-angle shadow evap- 
oration of Al and Cu through a polymer resist mask. 
First, an Al layer of thickness d$ ~ 30 nm was deposited, 
followed by in situ oxidation in the e-beam evaporator 
chamber in a few millibars of pure oxygen for a few min- 
utes. Finally, a Cu layer of thickness — 30 nm was 
evaporated at a different angle, forming the N island with 
four tunnel contacts to the Al electrodes. In addition, Cu 
replicas of the Al leads form large area tunnel junctions 
by partly covering the Al layer, serving as QP traps, al- 
beit of suboptimal performance [19] . 

The island electron temperature Tn is measured by 
biasing the probe junctions by a constant current 2th 5 and 
measuring the voltage drop Vth, calibrated against To at 
V = 0. The solid lines in Fig. [2] (a) show the measured 
Tn as a function of V at bath temperatures To between 
0.1 and 0.5 K, in zero field (red) and at B± = 3 G (blue). 
In the following, the minimum Tn along each curve at a 
fixed To is denoted by TN, m in, and the corresponding bias 
voltage by V^ pt . The strong influence of small B± on the 
cooling is evident: At eV ~ 2 A the maximum cooling 
ST = TMmin — Tn,o at each To increases by several tens 
of percents. The cooler QP current displays analogous 
behavior, Fig. [2](b). At the same time the optimum bias 
voltage V op t increases, Fig. [5](d), while heating at V > 




FIG. 2. (color online) (a) Temperature Tn of the N island and 
(b) cooler SINIS IV characteristic at several bath tempera- 
tures To as functions of the cooler bias voltage V, in zero field 
(red) and B± = 3 G (blue), (c) Relative minimum tempera- 
ture in zero field (red symbols) and at B± ~ 3 G (blue sym- 
bols). The solid and dashed-dotted lines show the calculated 
temperature reduction for various degrees of thermalization of 
the QPs and of the island phonons (see text), (d) Optimum 
bias voltage V op t vs B± at T = 285 mK, with the correspond- 
ing temperature drop in Fig. [I] (a), (e) Phonon temperatures 
in zero (red) and optimum (blue) field, required to reproduce 
the observed TN, m in- Dashed line shows T p h = Tq. 



V^pt diminishes. The cooling enhancement is symmetric 
in the applied field. 

The improved refrigeration is summarized in Fig. [2] (c) , 
where the symbols show the T -dependent relative mini- 
mum temperature in zero field and close to optimum B± . 
We observed the improved cooling also with Ag as the 
normal metal, in single NIS junctions with various grad- 
ually widening lead geometries close to the junction, and 
in a parallel field. In the latter case, the required fields 
were larger by an order of magnitude and dependent on 
the field orientation in the sample plane. 

The thin-film Al leads behave as a type II superconduc- 
tor, so that B± penetrates in the form of vortices [29] l30]. 
Based on a typical normal state resistivity p = 3.5 pQcm 
of our Al at 4.2 K [31], the elastic mean free path is 
I ~ 8 nm. With the Bardeen-Cooper-Schrieffer (BCS) 
coherence length £o — 1600 nm and the London pen- 
etration depth Al — 16 nm for bulk pure Al at low- 
temperatures, one obtains £ = 0.855(£o0 1 ^ 2 — 100 nm 
and A = A L (£ /0 1/2 - 230 nm for our A1 films witn tne 
Ginzburg-Landau parameter ft = A/£~2.4> 1/V2 and 
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the lower critical field for the bulk material H c \ < 100 G. 
As sketched in Fig. [I] (a), the S leads of the cooler junc- 
tions have an initial width of approximately 1 /am. At 
a distance of 1 /im away from the island, they widen 
to 2.5 /im width and continue for 15 jam before again 
widening to 10 /im width and connecting to large-area 
bonding pads further 350 /im away. The magnetic field 
below which vortices are completely expelled from a 
long and narrow S lead of width W is of the order of 
A) = $o/W 2 [32], where $ = h/(2e) - 2x lO" 15 Wb 
is the magnetic flux quantum. For a strip of width 
W = 10 /im, B ~ 0.2 G, whereas W = 2.5 /im results in 
Bo ~ 3.3 G. Taking into account the demagnetizing fac- 
tor l — n z ~ (0.5... 2) x 10 -2 of our films we can conclude 
that the initial increase in \ST\ in Fig. [I] (a) observed 
below 1 G and the turn-back that starts close to 3 G 
are consistent with vortex penetration into the wide and 
narrow parts of the lead, respectively. 

In Ref. 3 with large area NIS junctions, the increased 
sub-gap conductance at small B± could be directly asso- 
ciated with the fraction of vortices in the junction area. 
In contrast, we do not observe an increase in the cooler 
junction sub-gap current in the small fields. The ther- 
mometer junctions with narrower S electrodes are not 
considerably affected in fields B± < 10 G even at bias 
voltages close to 2 A, indicating that vortices exist only 
further away from these junctions. 

We assume quasiequilibrium with a local electronic 
temperature different from the bath temperature [14]. 
The cooling power depends on the temperatures Tn and 
Ts of N and S electrodes near the interface. They are 
found from the equation of heat balance in the N island 

2Q N (V, T N , T s ) = E N V N (T p 5 h - T N ) + P , (1) 

and of the heat conduction in each superconducting lead, 

V-[KVT(x)] = q(x), (2) 

with the boundary conditions ^fnAT'ffi) = Qs? and 
T(x) — » To at x — » ±00 for the right or left lead (A is a 
wire cross section). The temperature of the lead at the 
interface is Ts = T(x = 0). The first term in the rhs of 
Eq. ([!]) describes the heat transferred to phonons in the 
normal island. The island volume and electron-phonon 
coupling constant are Vn = -4n^n and En, respectively. 
P ^ 1 fW in Eq. (§ is a constant residual power due to 
imperfect RF-filtering of the measurement. The heat Qn 
extracted from the island through a single NIS junction 
and the heat Qs injected into an S lead by tunneling are 

Qn,s = — |r- I n s (E s )E N ,s [fr N (E N ) - f Ts (E s )] dE. 

(3) 

Here, E N = E - eV/2, E s = E, /t n ,t s (£) = 
l/[exp(i£/fcB^N,s) + 1] are the Fermi occupation factors, 
and ns(E) = E /\fE 2 — A 2 ]| is the normalized BCS den- 
sity of states (DOS). 



The rhs of Eq. ([2| is the power transferred from the 
unit volume of the superconductor into the (unbiased) 
normal trap with temperature To. Similarly to Eq. (J3|, 

= -2^T I ME)E[f T (E) - f To (E)} dE 
= [£(T)-£(T )]/r tT , 

where r^ 1 = l/(2e 2 JV ptr^s) is the time of relaxation to 
the trap, p tr being the trap/superconductor tunnel re- 
sistance of unit contact area and JVo denotes the normal 
state DOS at the Fermi energy per one spin projection, 
while £(T) is internal energy of the superconductor with 
the gap A at temperature T(x) [33]. In Eq. Q we as- 
sume that the electronic subsystem releases heat to the 
normal trap rather than directly to the phonon bath. In- 
deed, the trap relaxation rates are r^ 1 ~ 10 6 ...10 T s _1 
for the contact resistances p tT ~ 3 . . . 0.3 kQ x (jam) 2 
while the electron-phonon relaxation rate in aluminium 
is t^ 1 < 3 x 10 5 5 _1 for the experimental temperatures. 
Equation Q is obtained by averaging it over inhomo- 
geneities in A, assumed to have low areal density and 
short scale compared to the inelastic relaxation length, 
i.e., the scale of T(x) variations. The thermal conduc- 
tivity k and the heat current into the trap q are spa- 
tially averaged quantities k = ks (1 — r) + t and 
q = qs (1 — r)+#N r for a superconductor having a normal 
fraction r proportional to B± which models the presence 
of vortices. The thermal conductivity [28], fts, and the 
heat current, q$, of a superconductor at k^T <C A are ex- 
ponentially suppressed relative to their normal-state val- 
ues ft N = C cr N T and q N = [(7r 2 A^ ^|/3)(T 2 -T 2 )]/r tr ac- 
cording to fts/^N = (6/7r 2 )(A//c B T) 2 e- A / feBT and q s = 
(2^A) 3 / 2 iVo4 /2 [V / Te- A / feBT -V^e- A//CBTo ]Ar tr . Here 
C = (7r 2 /3)(/c B /e) 2 , and <jn is the Al normal state elec- 
trical conductivity. 

Because of exponentially small and qs at k^T <C A, 
the temperature Ts can be very sensitive to the vor- 
tex fraction r. For reference, the black dotted line in 
Fig. [2] (c) shows the calculated TN, m in as a function of 
To in the limit of perfect thermalization, Ts = T p h = To. 
The red and blue solid lines are the results of Eqs. ([T]) 
and ([2| for equilibrium phonons in the island, T p h = To. 
The red line is obtained with r = 0, so that the S elec- 
trode overheats at most all the way to the large-area 
bonding pad. Especially towards the lowest bath tem- 
peratures, the observed cooling in zero field is consider- 
ably weaker than the prediction of the model for perfect 
thermalization of S electrodes. To estimate the effect in 
the optimum field, we set r = 1 in the 10 jam wide elec- 
trode section, but keep r = in the narrower section. 
The result is shown as the solid blue line in Fig. [2] (c). 
In our samples, a considerable fraction of the N island 
is located on top of the S electrodes such that T p h can 
be essentially higher than Tq. As a worst-case estimate, 
we assume the N island phonons to be overheated to 
^ph = Ts. The result for the above two cases of full and 
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partial S electrode overheating are shown as the red and 
blue dash-dotted lines, respectively. The predicted influ- 
ence of the field is now stronger, and the measured cool- 
ing in both zero and optimum field is bracketed between 
the solid and the dash-dotted line. Approximate phonon 
temperatures T p h at V op t that reproduce exactly the ob- 
served TN,min are displayed in Fig. [2] (e). In both zero 
field (red, upper curve) and in the optimum field (lower 
curve), Tph reflects the temperature T$ of the QPs close 
to the junction, ranging from close to 0.8Ts at Tq = 0.1 K 
to around 0.95Ts at To = 0.5 K. We stress that the large 
field-induced improvement evident in Figs, [I] and [2] is ob- 
served because of the considerable S electrode overheat- 
ing in zero field. It causes also the significant increase 
in V^pt as B± is increased from zero: V^ pt is close to the 
ideal value 2(A — 0.66kBT^)/e [14] only at optimum B±. 



for the field- improved QP relaxation in the S leads. The 
work can provide means for optimizing the performance 
of superconducting nanostructures, and sheds additional 
light on the unsolved problem of nonequilibrium and 
residual quasiparticles. 
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FIG. 3. (color online) Cooling curves in selected magnetic 
fields for parallel SINIS coolers with different S electrode ge- 
ometries, at fixed To c± 130 mK. (a) When initially wide S 
leads are followed by a narrower section, applying a finite B± 
weakens the cooling monotonously, (b) In a structure with 
initially narrow leads, the behavior is non-monotonous, with 
enhanced cooling in small fields. The insets sketch one half 
of the cooler structure with 10 NIS junctions in parallel. 



To emphasize the role of the S electrode geometry, we 
performed additional experiments on parallel SINIS cool- 
ers. As shown in Fig. [3] (a), in a sample with initially wide 
leads, vortices form first close to the junctions, and apply- 
ing B± monotonously weakens the cooling. In contrast, 
with narrow leads close to the junctions as in Fig. [3] (b), 
the cooling is optimized at a finite B±. 

In conclusion, we observed that a small magnetic field 
enhances relaxation processes in a superconductor thus 
stabilizing superconductivity in non-equilibrium condi- 
tions. Significantly improved electronic cooling in a tun- 
nel junction was achieved in a small perpendicular mag- 
netic field. The enhancement of relaxation can be rele- 
vant also for "re-entrant superconductivity" observed in 
Zn nanowires [lOl [IT] driven out of equilibrium by su- 
percritical current. A quasiequilibrium model accounts 
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